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Abstract 

Many sea level rise assessments focus on populations presently inhabiting vulnerable 
coastal communities, but to date no studies have attempted to model the destinations of 
these potentially displaced persons. With millions of potential future migrants in heavily 
populated coastal communities, sea level rise scholarship focusing solely on coastal com¬ 
munities characterizes sea level rise as primarily a coastal issue, obscuring the potential 
impacts in landlocked communities created by sea level rise induced displacement. Here I 
address this issue by merging projected populations at-risk of sea level rise with migration 
systems simulations to project future destinations of sea level rise migrants in the United 
States (U.S.). I find that unmitigated sea level rise is expected to reshape the U.S. pop¬ 
ulation distribution, potentially stressing landlocked areas unprepared to accommodate 
this wave of coastal migrants - even after accounting for potential adaptation. These 
results provide the first glimpse of how climate change will reshape future population dis¬ 
tributions and establishes a new foundation for modelling potential migration destinations 
from climate stressors in an era of global environmental change. 
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1. Introduction 

It is generally understood that sea level rise of lrn to 2m (DeConto and Pollard 2016; 
Vermeer and Rahmstorf 2009; Levermann et al. 2013) could lead to widespread human 
migration (IPCC 2014; Strauss, Kulp, and Levermann 2015) as residents of highly vulnerable 
coastal communities look to escape rising water levels. With up to 180 million people directly 
at risk to sea level rise in the world and over 1 billion living in the lower elevation coastal zone 
(Nicholls et al. 2011; Neumann, Vafeidis, et al. 2015), understanding the ramifications of 
these potential migrants on destination communities is a priority for climate change research 
(Black, Bennett, et al. 2011; Doos 1997; Mueller, Gray, and Kosec 2014; Hugo 2011a). 

Sea level rise assessments, identifying both the number and locations of potentially 
displaced persons, fill the literature (Strauss, Kulp, and Levermann 2015; Hauer, Evans, and 
Mishra 2016; Ben and Anders 2014) and are useful for the deployment of critical infrastructure 
in coastal areas. Yet questions of where the millions of potentially displaced persons will go 
remain unanswered despite a general understanding that sea level rise displaced persons are 
likely to have profound effects on future population landscapes (Doos 1997; Ben and Anders 
2014). Only a few studies have put forth general hypotheses regarding sea level rise migration 
(Curtis and Schneider 2011; Doos 1997) and this void has prompted recent calls for additional 
migration modeling (Willekens et al. 2016; Hugo 2008). To date, no studies modeling precisely 
how sea level rise induced migration will affect the population distribution exist. By focusing 
solely on coastal communities without directly addressing sea level rise induced migration, we 
likely underestimate the scale, magnitude, and geographies of these impacts. 

Relationships between environmental stressors and migration are highly complex as 
press and pulse events trigger migration responses that range from short-distance temporary 
migration to permanent long-distance migration; some will move and others will not (McLe- 
man 2013a; Gutmann and Field 2010; Collins et al. 2010; Hunter, Murray, and Riosmena 
2013; Thiede and Brown 2013; Entwisle et al. 2016). Sea level rise is unique among environ¬ 
mental stressors as the conversion of habitable land to uninhabitable water is expected to lead 
to widespread human migration without the deployment of costly protective infrastructure 
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(Doos 1997; IPCC 2014; Curtis and Schneider 2011; Strauss, Kulp, and Levermann 2015). It 
is unclear, however, what will actually trigger future climate migrants: press events, such as 
drought or sea level rise, etc. or pulse events such as tropical cyclones. When climate effects 
are integrated over long periods of time, it is likely that a combination of press and pulse 
events will spur migration (Arenstam Gibbons and Nicholls 2006) across preexisting migra¬ 
tion pathways (Hunter, Murray, and Riosmena 2013; Gutmann and Field 2010), leveraging 
established networks of social capital and kin networks in destination decisions (Findlay 2011). 
This is because press and pulse events that spur migration operate mostly independently of 
the kin networks and social capital that drive destination decisions (Black, Bennett, et al. 
2011). Thus, climate migrants resulting from press stressors will likely constitute “enhanced”, 
or extra, normal out-migration. 

I combine estimates of the populations at-risk to sea level rise (Hauer, Evans, and Mishra 
2016) within a migration systems simulation to estimate both the number and destinations 
of potential sea level rise migrants in the United States (U.S.) over the coming century. By 
focusing on the destinations of sea level rise migrants I am able to more holistically describe 
the impacts of sea level rise. This study aims to answer one fundamental question regarding 
sea level rise induced migration: Where will sea level rise migrants likely migrate? Local 
officials in landlocked communities can use these results to plan for potential infrastructure 
required to accommodate an influx of coastal migrants and could shift the conceptualization 
of sea level rise from a coastal issue to an everywhere issue. 

2. Methods and Materials 

I utilize two primary sources of data concerning the magnitude of flows in the migration 
system and the migration system itself. The first, the magnitude of flows, comes from pub¬ 
lished populations projected to be at risk to sea level rise (Hauer, Evans, and Mishra 2016). 
Hauer et al. (Hauer, Evans, and Mishra 2016) used the National Oceanic and Atmospheric 
Administration’s (NOAA) 0m, 0.9m (3 feet), and 1.8m (6 feet) sea level rise datasets for 
twenty-two coastal states and the District of Columbia. These datasets simulate expected 
changes in the mean higher high water (MHHW) mark on areas that are hydrologically con- 
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nected to coastal areas without taking into account additional land loss caused by other nat¬ 
ural factors such as erosion. They projected populations using a modified Hammer Method 
(Hammer et al. 2004) combined with the Housing Unit Method (Hauer, Evans, and Alexander 
2015) for population estimation to create temporally contiguous sub-county boundaries over 
a 70-year base period from 1940-2010, which were then used to project populations at these 
same sub-county geographies through the use of a linear/exponential extrapolation approach 
for projecting Census Block Groups (CBG) from 2010 to 2100. The populations at risk to sea 
level rise, aggregated to the US County and available in their Supplementary Table 2 (Hauer, 
Evans, and Mishra 2016), provide the magnitude of out migrants from 319 coastal counties. 

For this research, I used the Internal Revenue Service’s (IRS), produced in conjunction 
with the US Census Bureau, annual series of county-to-county migration datasets (IRS 2016) 
as the basis for the migration system. The IRS datasets utilize the IRS Individual Master 
File containing every Form 1040, 1040A, and 1040EZ processed by the IRS and includes 95 
to 98 percent of all individual tax filers and their dependents. The Census Bureau identifies 
migrants when a current filing years’ return is from a different location than the matched pre¬ 
ceding years’ return. These data only capture the tax-filing universe but the spatio-temporal 
stability (Rogers, Little, and Rayrner 2010; Engels and Healy 1981; Franklin and Plane 2006) 
coupled with the very large administrative sample make them attractive for modeling large 
scale migration patterns. The IRS does suppress migration flows comprising fewer than 10 
individual migrants, systemically suppressing small rural migration flows. However the long¬ 
term trend of rural out-migration to urban areas (Lichter and Brown 2011) is expected to 
continue in this century (Jones and O’Neill 2013). 

Migration systems theory (MST) has been tied to environmental migration in recent 
years (DeWaard, Kim, and Rayrner 2012; Fawcett 1989; Curtis, Fussell, and DeWaard 2013; 
Fussell, Curtis, and DeWaard 2014). MST is a branch of migration research that holistically 
examines migration options by studying all origin-destination combinations rather than any 
single origin-destination combination (Fawcett 1989; Massey et al. 1994; DeWaard, Kim, and 
Rayrner 2012). Migration decisions - not just the decision to migrate, but also decisions on 
where to migrate - are often driven by kin networks, employment opportunities, amenities, 
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Figure 1: Sample migration systems from the IRS’ county-to-county migration 
dataset for Eagle County Co, Assumption Parish LA, and Dougherty County 

GA. The thicker the band, the more migrants. Blue bands represent destination systems 
while red bands represent origin systems. 


both natural and economic, economic vitality, and so forth (Lee 1966; Pandit 1997; Fawcett 
1989; Haug 2008; Thiede and Brown 2013; Black, Bennett, et al. 2011). This network of “pull” 
factors embedded within the migration system tends to drive locational decision making due 
to environmental, or other, “push” factors (Hugo 2008; Hugo 2011b; Schultz and Elliott 2013; 
Gray and Bilsborrow 2013; McLeman 2013b; Findlay 2011; Black, Bennett, et al. 2011). 

To describe the complete migration system in the United States, let matrix M{x ) rep¬ 
resent all possible county level origin-destination combinations. The sum of any given column 
and row in the matrix will equal the total number of migrants into or out of any given county. 
For this analysis I am only concerned with the 319 coastal counties expected to experience 
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some form of sea level rise inundation under the 1.8m scenario and their connections to the 
other 3,113 US counties (n=993,047 matrix cells). I created these matrices for each year of 
the IRS migration data between 1990 and 2013. Figure 1 shows examples of these systems 
into and out of three sample counties. 


M(x) = < 


M i4 


M. 


3,i 


where 


o € 1, ...,319 
d € 1, ...,3113 


Mi, 


M- 


= K,d] 


3,3 


I employed the use of an unobserved components model (UCM) for forecasting equally 
spaced univariate time series data (Harvey 1990). UCMs decompose a time series into com¬ 
ponents such as trends, seasons, cycles, and regression effects and are designed to capture 
the features of the series that explain and predict its behavior. UCMs are similar to dynamic 
models in Bayesian time series forecasting (Harrison and West 1999). All projections were 
undertaken in SAS 9.4 using the PROC UCM procedure. 

The basic structural model (BSM) is the sum of its stochastic components. Here I use 
a trend component /at and a random error component e t and it can be described as: 


Ut = 1H + et 


( 1 ) 


Each of the model components are modeled separately with the random error St modeled 
as a sequence of independent, identically distributed zero-mean Gaussian random variables. 
The trend component is modeled using the following equations: 
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IH = IH- 1 + Pt-i + Vt 
Pt = Pt -1 + 

Vt ~N(Q, <tJ) 

6 ~ iV(0, of) 

These equations specify a trend where the level and the slope Pt vary over time, 
governed by the variance of the disturbance terms r/t and £< in their equations. Here all 
origin-destination dyadic pairs containing any migration information over the series were 
modelled (n=46,203) in individual UCM models. 

This approach allows for the projected evolution of the migration system as dyadic 
pairs either strengthen or weaken over time allowing for migration links to wax or wane over 
the projection horizon. Empirical simulations of environmental migration have proven very 
fruitful in the modeling of climate induced migration (Reuveny and Moore 2009; Curtis and 
Schneider 2011; Feng, Krueger, and Oppenheimer 2010) and here I build upon those efforts 
by projecting future climate migration. 

Our current understanding of the migratory response to sea level rise is still underde¬ 
veloped. Will displaced coastal populations relocate into the parts of coastal communities 
unaffected by sea level rise? Or will displaced persons migrate to more inland areas free 
from the challenges of sea level rise (Doos 1997; Curtis and Schneider 2011)? Many areas in 
threatened coastal communities are still eligible for human settlement and could be possible 
destinations for future sea level rise migrants. To capture both possibilities, I employ a rak¬ 
ing procedure to proportionally adjust in-migrants based on the inverse of the proportion of 
the population affected by sea level rise and a redistribution of those migrants to unaffected 
counties. 


Mid = Ml d • (1 - i^/Pj) 


(2) 
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The adjusted number of in-migrants M to destination county d from origin county o 
at time t is equal to the number of migrants multiplied by one minus the proportion of the 
population impacted by sea level rise in the destination county at time t. For inland counties, 
the right side of the equation equals one, yielding no adjustment. 


Mod{i=U] 


(E m »-E*o) 


/ Mpd[i=U] \ 


( 3 ) 


However, the out-migrants from each affected coastal county must be raked to equal 
the total population at-risk to sea level rise expected to be displaced. This is accomplished by 
redistributing the difference of the unadjusted migrants M f 0 from origin o from the adjusted 
migrants from origin o and multiplying it by the proportion of adjusted migrants from origin 
o to destination d to the unaffected counties (i = JJ) from the total adjusted migrants from 
origin o. In this way, the underlying migration system from each individual origin is preserved 
in the raking procedure. 

It is likely that many communities will deploy a wide variety of adaptation measures 
including sea walls, beach and marsh nourishment, pumps, or elevate homes and roads to 
protect both people and property and IPCC reports have increasingly emphasized adaptation 
when discussing sea level rise (Brown et al. 2014). Global estimates of adaptive infrastructure 
for sea level rise could reach US$421 billion (2014 values) per year (Nicholls et al. 2011) and 
could cost upwards of US$1.1 trillion in the U.S. (Neumann, Emanuel, et al. 2015). However, 
the deployment of adaptation measures is driven by wealth for both cities (Georgeson et al. 
2016) and individuals (Black, Bennett, et al. 2011; Adger, Arnell, and Tompkins 2005). To 
approximate this dynamic, I assume that households earning greater than US$100,000 per 
year in the 2011-2015 American Community Survey are likely to adapt to sea level rise in 
some manner and thus unlikely to migrate. This income threshold represents approximately 
the top quartile and double the U.S. median household income and is neither too restrictive 
nor too broad to capture the range of individual adaptive measures. 
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2.1. Evaluation of Migration System Projection 

Table 1: Number of IRS migration counts that fall within the 2/3 projection 
interval. This table contains the percentage of dyadic origin-destination pairs where the 
projected migration flows fall within the 2/3 projection interval. 


Year 

# within the 2/3 

Projection Interval 

# of Dyadic 

Migration Pairs 

Percent 

2004 

34,719 

46,203 

75.1% 

2005 

32,391 

46,203 

70.1% 

2006 

33,458 

46,203 

72.4% 

2007 

33,566 

46,203 

72.6% 

2008 

33,960 

46,203 

73.5% 

2009 

34,828 

46,203 

75.4% 

2010 

34,530 

46,203 

74.7% 

2011 

32,217 

46,203 

69.7% 

2012 

31,836 

46,203 

68.9% 

2013 

33,075 

46,203 

71.6% 

TOTAL 

334,580 

462,030 

72.4% 


Projection intervals allow us to examine the feasibility of the future projected migration 
systems and are typically employed in the evaluation of population projections (Gerland et 
al. 2014). Population geographers have typically used the 2/3 or 66% projection interval 
to assess the accuracy of a population projection (Tayrnan, Smith, and Lin 2007; Swanson 
and Beck 1994), representing “low” and “high” scenarios that are ’’neither so wide as to be 
meaningless nor too narrow to be overly-restrictive.” (Swanson and Tayrnan 2015). 

To examine the feasibility of the migration system projections, I produce projections 
based on the equations in the preceding section with the base period 1990-2003 and an eval¬ 
uation period of 2004-2013. If less than 2/3 of the IRS migration counts fall within the 2/3 
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projection interval then the results would suggest less than ideal accuracy. However, if more 
than 2/3 of the IRS migration counts fall within the 2/3 projection interval, it would suggest 
an ideal amount of accuracy. I assessed the 2/3 interval for 10 years of projections for an 
evaluation of base period 1990-2003 and a projection period of 2004-2013. 

Table 1 shows the overall number of IRS migration counts that fall within the 2/3 
projection interval for each year of the evaluation period. Overall, the UCMs produce robust 
projections as all projection years are above the 2/3 projection interval. 


3. Results 

I find that in the U.S., every state, 86% of U.S. Core Based Statistical Areas (CBSA) 
(791 out of 915), and 56% of counties (1,735 out of 3,113) could be affected in some way 
by net migration (in-migration minus out-migration) associated with 1.8m of sea level rise 
(Figure 4). Florida could lose more than 2.5 million residents due to 1.8m of sea level rise 
while Texas could see nearly 1.5 million additional residents (Figure 3). Additionally, nine 
states could see net-losses in their populations due to sea level rise. Figure 2 demonstrates 
all origin-destination flows at the U.S. state level demonstrating that the sheer magnitude of 
places affected could alter the U.S. population landscape, sea level rise migrants are expected 
to comprise both intra- and inter-state migrations and no state is left untouched by sea level 
rise migration. Even accounting for the deployment of adaptive infrastructure, millions of 
people could still migrate (Table 2). 

My results also suggest that CBSAs such as Austin TX, Orlando FL, Atlanta GA, 
and Houston TX could see more than 250,000 previously unforeseen future sea level rise 
net migrants each (Table 2). Thirteen CBSAs could see more than 100,000 sea level rise 
net migrants by 2100 with 1.8m of sea level rise. Conversely, ten CBSAs could lose more 
than 100,000 residents due to sea level rise with Miami FL losing over 2.5 million residents. 
Even accounting for those who could adapt in place, many inland communities could see 
tens of thousands of sea level rise in migrants and many coastal communities could lose tens 
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Figure 2: Circular plot of bi-lateral sea level rise migration flows for US States 
under the 1.8m scenario and no adaptation. Tick marks show the number of migrants 
(inflows and outflows) in thousands. States are ordered clockwise by the size of inflows. 
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Table 2: Select core based statistical area (CBSA) destinations of sea level rise net 
migration (in-migration minus out-migration) in 2100 with and without adapta¬ 
tion under the 1.8m scenario. +/- represents the 80th confidence interval. 


CBSA 

No Adaptation 


With Adaptation 


Net-Migration 

+/- 

Rank 

Net-Migration 

+/- 

Rank 

Austin-Round Rock, TX 

818,938 

243,821 

i 

625,627 

179,186 

i 

Orlando-Kissimmee-Sanford, FL 

461,411 

62,665 

2 

369,120 

38,834 

2 

Atlanta-Sandy Springs-Roswell, GA 

320,937 

131,984 

3 

248,684 

68,868 

3 

Phoenix-Mesa-Scottsdale, AZ 

100,524 

12,851 

13 

73,935 

1,949 

13 

Myrtle Bcach-Conway-N Myrtle Beach, SC-NC 

12,146 

13,855 

78 

3,142 

9,389 

141 

North Port-Sarasota-Bradenton, FL 

-208 

25,057 

727 

2,128 

21,717 

168 

Los Angeles-Long Beach-Anaheim, CA 

-3,140 

51,590 

737 

13,181 

22,200 

61 

New York-Newark-Jersey City, NY-NJ-PA 

-50,804 

494,625 

775 

15,808 

194,047 

50 

New Orleans-Metairie, LA 

-500,011 

24,053 

795 

-373,283 

10,733 

795 

Miami-Fort Lauderdale-West Palm Beach, FL 

-2,509,978 

155,119 

796 

-2,009,263 

95,845 

796 


of thousands of residents. Extended results for all counties and CBSAs are available in the 
supplementary materials (Supplementary Dataset). 


With many projected migrants remaining in coastal communities (Table 2 and Sup¬ 
plementary Dataset), sea level rise could generate millions of “trapped” people (Black et al. 
2013). Trapped populations are sometimes discussed through the concept of involuntary im¬ 
mobility (Black et al. 2013), but there is also those who do not desire to move and thus 
constitute voluntary immobility. These results suggest that many people displaced by sea 
level rise could find themselves or their descendants exposed to sea level rise, even with mi¬ 
gration as an adaptation, as sea levels continue to rise past the year 2100 and migration that 
constitutes relocation to presently safe, but ultimately vulnerable, coastal communities. 


Additionally, infrastructure challenges required to protect coastal communities are well 
documented (Hauer, Evans, and Mishra 2016; Strauss, Kulp, and Levermann 2015), but the 
infrastructure challenges of accommodating millions of sea level rise migrants in largely unpre¬ 
pared inland municipalities is virtually unexplored. For many destinations, such as Riverside 
CA, Phoenix AZ, Las Vegas NV, and Atlanta GA, already experiencing water management 
and growth management challenges, the sea level rise migrants who wash across the landscape 
over the coming century could place undue burden in these places if accommodation strate¬ 
gies are left unplanned. Studies of migration impacts do not solve the challenges in these 
areas, but rather reveal a more holistic understanding of sea level rise impacts and needed 


interventions. 
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Figure 3: Net change in population due to sea level rise under the 1.8m scenario 
and no adaptation. I considered migration destinations for all 50 states and the District 
of Columbia and migration origins for 22 states and the District of Columbia. These are the 
net changes in population due to both in- and out-migration due to sea level rise. 
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Figure 4: Estimated sea level rise net migrants (in-migrants minus out-migrants) 
for counties (a) and core based statistical areas (b) under the 1.8m scenario and 
no adaptation. I considered migration destinations for all 50 states and the District of 
Columbia and migration origins for 22 states and the District of Columbia. These are the net 
changes in population due to both in- and out-migration due to sea level rise. 

4. Discussion 

sea level rise has been broadly conceptualized as a coastal issue or hazard, as assess¬ 
ments have focused on the effects in coastal communities (Strauss, Kulp, and Levermann 
2015; Hauer, Evans, and Mishra 2016; Neumann, Vafeidis, et al. 2015; Nicholls et al. 2011). 
With millions of potential future migrants in heavily populated coastal communities, sea level 
rise scholarship focusing solely on coastal communities endorses a narrative that character¬ 
izes sea level rise as primarily a coastal issue, obscuring the potential impacts in landlocked 
communities created by sea level rise driven migration. My work shows that this coastal 
conceptualization of sea level rise creates a deceptively small area of affect if relocation is left 
unaccounted. This work offers the first glimpse of how sea level rise could alter the popu¬ 
lation distribution of the U.S. as both coastal and landlocked communities are likely to be 
affected by sea level rise: directly in coastal areas due to sea level rise itself and indirectly in 
landlocked areas through the influx of people escaping sea level rise. 

Furthermore, the migration approach for examining destinations associated with cli¬ 
mate change shown here allows for modeling migration destinations of other climate change 
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stressors. For instance, it has been estimated that parts of the Middle East and North Africa 
(MENA) could become uninhabitable by the end of the century, potentially spurring an exo¬ 
dus of 500 million people (Lelieveld et al. 2016). Future scholars could employ my approach 
to model the destinations of these potential MENA migrants. There is tremendous poten¬ 
tial in coupling migration systems information with climate change models to examine the 
implications of climate change induced migration. This type of modeling requires detailed 
origin-destination migration information, limiting the areas where this approach could be used 
to those where data are available. 

Migration models are only as good as the data underlying them. The suppression of 
IRS migration data of flows with fewer than 10 migrants could systematically bias my results 
against rural areas far from coastal communities and the modeling approach undertaken does 
not allow for new unforeseen migration pairs to emerge in the future. However, the destination 
counties cover 93.6% of the U.S. population and over 250 destination counties are outside of 
CBS As (the typical distinction of urban/rural), limiting the scope of the rural bias to mainly 
sparsely populated communities far from coastal areas. New destinations for sea level rise 
migrants could still emerge in rural areas, and if they do, my results of the geographic spread 
of sea level rise migrants could be considered conservative. Environmental migration scholars 
could further investigate the likelihood of emerging rural destinations related to climate change 
migration to better inform future modeling efforts. 

Previous trends are not always indicative of future results. Societal and economic shifts, 
population ceilings, local growth ordinances, adaptive behavior, and climate change itself 
could all change future migrant destinations. Although I model the potential destinations of 
sea level rise migrants, I do not precisely model how other climate stressors or other factors 
might influence the future migration systems. Our current understanding of the location 
decisions of environmental migrants is still limited and there have been recent calls to better 
understand migration flows (Willekens et al. 2016). My approach builds on the growing 
literature concerning migrant destinations and environmental change (Findlay 2011; Smith 
et al. 2011; Black, Kniveton, and Schmidt-Verkerk 2011; Black, Bennett, et al. 2011; Black, 
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Adger, et al. 2011; Curtis and Schneider 2011; Doos 1997) and accounts for potential future 
migration systems. 

High potential exists for development of deeper and more integrated examinations of 
the role of adaptation and other climate stressors on SLR migration. Neglecting to account 
for adaptive behavior could lead to an overestimation of either flood risk (Haer et al. 2016) or 
migration (Kocornik-Mina et al. 2016). Future studies could examine varying adaptation sce¬ 
narios related to the 100-year flood plain or the lower-elevation coastal zone creating differing 
scenarios of “stayers.” 

If future migration pathways mimic past pathways, sea level rise is expected to reshape 
the U.S. population distribution and could stress some landlocked areas unprepared for these 
migrations while revitalizing others, sea level rise is currently framed as a coastal hazard but 
the migratory effects could ripple far inland. My results show the importance of accounting 
for future migrations associated with climate change in long range planning processes for 
disaster management, transportation infrastructure, land use decisions, etc.. 

Data Availability 

The data that support the findings of this study have been deposited in openICPSR 


(10.3886/E100413V3). 
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